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ABSTRACT 

A strong constraint on 13th order (odd degree) terms in the geopotential has 
been derived from Navy tracking on a DIADEME 2 fragment (1967-14F). This 
object (perigee height: 580km, orbit inclination: 38. 9“) is presently decaying 
slowly through perfect commensurability with these terms. The resonance 
forces will increase its inclination by 0. 02° when the passage is complete by 
late 1974. The constraint (lumped harmonics), derived by adjustment of a pair 
of harmonic coefficients to the Navy inclination data (principally) is: 

10 9 (14.8 ± 0.8, 48.3 ± 0.7) = 0. 023(C, S) 1313 - 0.172(C,S) 1513 + 0. 505(C,S) 17>13 
-0.884(C,S) 19>13 + (C,S) 21>]3 - 0.673(C,S) 23>1 3 + 0.099(C,S) 25 ,i 3 + 0.295(C,S) 27 ,i3 

-0.279(C,S) 29tl3 + 0.018(C,S) 3U3 + . . . 

There should be a significant contribution to this result from terms as high as 
29th degree. But current geopotential solutions (for 13th order terms) to this 
degree are about 20% in error when judged by this independent data. 
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13TH ORDER RESONANCE FROM NAVY TRACKING ON A 


DIADEME 2 FRAGMENT 


INTRODUCTION 

Among the thousands of Earth satellites launched in the past 15 years, hun¬ 
dreds must have suffered or will suffer measurable effects from poorly known 
resonances with the Earth’s geopotential (King-Hele, 1973a; Gabbard and 
Wackernagel, 1971). But less than 50 of these objects have actually been used 
geodetically. The major difficulty of course is the lack of adequate tracking for 
the majority of the objects which are (or were) rocket bodies or fragments 
(debris) of larger satellites. However, in the case of deep resonance or close 
commensurability to the Earth’s rotation, the effects are so large that even 
tracking with crude instruments can reveal them (King-Hele, 1973b). Many of 
the observations made in support of the tracking on decaying 15th order resonant 
orbits were accomplished with amateur observers using binoculars and 
stopwatches. 

In these deep resonances the change of the inclination of the orbit is of the 
order of 0. 01°. To achieve overall inclination accuracies of a tenth of this fig¬ 
ure, a topocentric angle good to only about a minute of arc is necessary. King- 
Hele's results were often much better than this, but his amateurs were backed 
up by more precise camera observations. (See also Winterbottom and King- 
Hele, 1974.) In addition, King-Hele’s group has utilized minute-of-arc data 
from the U.S. Navy skin track Naval Space Surveillance (NAVSPASUR) system, 
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of which there is a great abundance on low altitude satellites. In fact since 
1972, the NAVSPASUR data and mean elements derived from them (by the Navy) 
dominate the orbits used in the 15th order resonance analysis (i.e., King-Hele 
and Walker, 1972). This data consists of (direction cosine) pointing angles 
from a fence of receiving stations across the United States. The signals are 
reflections from the satellite of a continuous high power radar fan illuminating 
space across the U.S. at about 28° latitude. The system has tracked objects as 
high as 18,000km altitude. Providing the inclination of the orbit is greater than 
about 28° and the object is large enough to be tracked, data on each revolution 
will be received. 

The existence of this good data on a large number of objects opens up the 
possibility of examining in detail all the deep geopotential resonances besides 
the one’s of 15th order. In this study, the Navy’s own mean elements were 
used directly. They proved sufficiently accurate to reveal the strong effect on 
the inclination of the decaying (13 revolutions/day) orbit of a DIADEME 2 
fragment (1967-14F). 

ANALYSIS 

In Figure 1, the Navy Mean Inclinations are shown from February 1972 to 
January 1974. Also in Figure 1 is a plot of the mean (primary) resonant longi¬ 
tude rate for this 13 revs/day orbit: 

where oj, M and flare the orbit’s argument of perigee, mean anomaly and right 

4 

ascension of the ascending node, and 8 is the Earth’s rotation rate. 
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Figure 1. U. S. Navy Orbit Inclinations for 1967-14F 



It is recalled that the longitude rates which determine the frequency of gravi¬ 
tational perturbations on an orbit in Kaula’s development of the geopotential 
(Kaula, 1966; especially p. 40, 49 and 55) are given as: 

p = ( £ - 2 p) to + (£ - 2 p + q) M + m (£l - 6 ), 
where £ and m are the degree and order of a geopotential harmonic, and p and q 
are additional indices related to the inclination and eccentricity of the orbit. 

The resonances (ip = 0, m 4- 0) will occur for (£ - 2p + q) M close to m0, 
since to and tl are small. The primary resonances (strongest) occur when # - 

2p + q = 1, since then m will be minimum, the closest integer to M in 

♦ 

revolutions/day. Secondary resonances (sub harmonics) for a given M (near a 
rational number of revs /day) will occur for £-2p + q = 2,3, 4. . . , but in 
any case all the resonant longitude rates can be characterized by the order m 
and the q index, by writing 

„ = -qto + (£- 2p + q) (to + M) + m {tl - 
For a given resonant order m, there will be a series of "side band” resonances 
characterized by q around the (generally) dominant one for q - 0. The q = 0 
resonance is also the mean of the series since q can take on all positive or 
negative integers. 

■ 

It is seen in Figure 1 that the longitude rate 13,0 for 1967-14F goes to 
zero over the period of record as the inclination of the orbit suffers a progres¬ 
sive oscillation of increasing period and amplitude. The increasing period 
closely matches the period of \I/ l3>0 . There is clearly a strong resonant 
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perturbation of this orbit entirely analogous to the decaying 15th order orbits 
first analyzed by Gooding (1971) and the 11th order orbit examined by Wagner 
(1973). 

But there is more significant detail in Figure 1 than resonance. The ''raw" 
mean inclinations given by the Navy (see also Table 1) are actually Brouwer ele¬ 
ments (determined over independent 7 day arcs) with both short and long period 
zonal effects removed (Brouwer, 1959). Only J 2 to J 5 are used in the Navy- 
Brouwer model and the coefficients are not up to date. The period of the princi¬ 
pal odd zonal inclination effect (2 tt/cj) is 65 days. The amplitude, with the cor¬ 
rect model, is about 0. 005°. The error in the Navy model is certainly seen in 
the raw inclination data. But, in addition, there is an 84 day lunar perturbation 
with amplitude 0. 0014° which is also observable. 

To clarify the quality of this "signal" I have added back the long period zonal 
effects used in the Navv-Brouwer model to produce mean elements free of this 
bias (also Table 1). Then I compared these (less smooth) mean elements to 
values calculated from a trajectory which includes all significant (but nonreso¬ 
nant) long period effects on the orbit, in particular the zonal perturbations from 
the Smithsonian (SAO) Standard Earth 2 (Gaposchkin and Lambeck. 1971). radia¬ 
tion pressure, atmospheric drag and direct lunar-solar gravity. The compari¬ 
son (observed minus calculated values) in Figure 2, shows the resonant signal 
much more strongly and also reveals what appears to be a much reduced resid¬ 
ual effect due to odd zonal error in the Smithsonian field. 
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C* 62072300000-01 
— y i j ! 704 4 6 42dO ~Q 1 


INCLINATION 

(DEGREES) 

I 

l' 

3B,92290000 
3£U 9323.93-3.4- 
38* 9214du0o 
^3325G-9-4- 


ARGUMENT OF 
PERIGEE 
(DEGREES) 

GJ 

f 

W 

17 1*4 370^ 00 
174^-404^ 0*2- 
210 *060 1 j 00 
31 0 -* 1 ' 9- ' ^-2.» !« * 6 4 


RIGHT ASCENSION OF 
THE ASCENDING 
NODE (DEGREES) 

£2 

&' 

2750 545 104)0 
. 2-7- 5^-9 4W34 5— 


MEAN ANOMALY 
(DEGREES) 

M 

M' 

190,6 1264000 
- 3 - fr Oa- S 58 6363- 


245«d750 GO 0 
e4 b 'o 86^660 4- 


14464607000 
-H4i r»i e e sa ? - 


3 * 1i?790210 
1^.137-79-021^- 
1 * 1£77c3 2y0 
i * 1277csZdO 
1 • 167705240 


0*62037700000-01 

C « U2001OOOOOD-U1 
C * 31-4 2-1 2337 
C * d1967300000-01 
.£-*.31-0 5 3* 3 60 60. - 01 - 


3d*9227J000 


36.91dluOOO 


24 9 *607 2j Oy 
..24 ^-«4ol> J - k - 0 -3- 


3d i 


91940000 
-92 1 323-UL- 


2b9 *072 1> 00 
-2o 9 l* 2 0-i 5>-3 9 — 
326 * 234 2J 0 0 
3-2 -8- . -P db -4-t-o -a — 


£ 15*6033000 
--21^60 C 7 - 5G1 - 
185.3944000 
-l-8S^3967-&99 


101 o45 17-0 00 
14 U 5924960- 


155 *3127000 
1-5- 5 * 3 1 6 9669- 


€2*28500000 
-62^16299356— 
27.12870000 
26 -u r7 7 6 1 d - 03 -3— 


l * Id77c413C 
1 * 167734 130 
1 *1 c 7 7 b 3 1 72 
1 * 137733170 
1 . 1 5776-J76C 
1.1877L3760 
1 * 187763330 
l .1 j77<Jib3‘0 
1 * 187731CcO 
1 o 1 c7 7 < 5 1 C 3 0 
t * 1C7779390 
1 *1£777>390 
1 ai37779330 
i♦1c7779350 
1 * 1 L 777 c330 
1 * 1677 7o35G 
1 « I c 7 77 5870 
1 *l£7 7?557C 
1 * l 67772430 
1 * 1277 72450 
L * 1=7765240 
1 * 157 7 6> 24 0 
1 . 1 t77fcc 490 
i * 147766*4 90. . 
1 * 167762530 
1«1 t 7762330 
1 * 1t7752 2 JO 
1«lc775423t 
1 * 16 77 34 24 C 
1 * 1 = 775424 0 
\ .1-7752 09 0 
1 *1L775209G 
1 ♦ 1 4774q 740 
.1* la? 7.4 6 74.0._ 


0* 51914600G0O-O1 
0 * dl 9-3-7553460-^01- 


36*9172JU00 
23xlo42L67JL3iJ- 


C*31750400000-01 
0.322-02 7.3.9.950 -.03 — 
G • 5 i7 8E5000OD— 01 
.0 * 82 331 7.56 330=01— 


38*91570000 

3a±**13J*2-±2-. 


38* 91580000 
:\ 8 A 91 1 9 3 9 3 9 


7 #33b?OOJ 00 
*7-663 9-4-. 3 2— 
46 *9320 0^ JO 
_4 7_*-2 dG-1.1^0 7- 
36*0701OJ00 
-<h6-*-J-U-j. 8 3 j<- 0 6— 


125*2260000 

125*2233766 


353 « 861000 0 
3-53o 4290419 


94*6 166 C 000 
94^6 217-2-09-0- 


2192.6943000 

-31^3652186 


64 * ?4l5CGQ0 
64 o 742 0 1 660 


203* 31 59000 

5 J - 


0»8l‘76060j000-01 
G* 322 67 2 772 7D-U.1-. 
C.61825200000-01 
G* dl98104 465D-01 
C* el 86 9 300 800—01 
-C . d 1 6 l 8 6d43.lD.-lU .. 
G.618131 00000-01 
G *31 2701622 30-01 
C.8166130000D-01 
C.61263593O00-01 
C. 8 16 4 48000OD-01 
QoOl 4oS.^3.475O-0L . 
C-dl757400000-01 
C*dl 7g03l3bdD — 01 
0*616938000 00-01 
Co-32 0 501 765q>0-ol 
C* 616802000 00-01 
.3)^32264219 flflaP-fl. 
C*616533 0000D*^8l 
G«822CG3c<+34D-01 
0* 61 60 360 OOOD — U1 
0*d19117812OD-03 . 
C * 61 626700000-01 
J fl 3143ol 92 72D-01 
C *bl 6 6u 78000D-01 
0 p.dl 1469 90.060-01._ 
C « ol 7 02o UGOOD-Ol 
C _«.8 1.0 9 3 4.0 0T_9D -C 1__ 


38*91550000 
.3 3^91 25 3 its ?_ 
38* 91510000 
3 6* 9141.66 9 0_ 
38*^1590000 
_3 0 *1 9.1 iLJLaf-b 1 j6— 


125* 30 57^ 0 0 
-12 3-*4) 43-4 *J94— 
164.8857u 00 

_16 4.*>+63 3^-27 — 

203 *99 55y 00 
—2 0 3* 615 3y-5-l—. 


34*65960000 

34-06-5-5-2873.5 


4 *258400000 


242*7033000 
242o9663684 
197*7912000 
193^2149702- 


2 34 * 18810CO 
3 34-0-16-136^3- 


152.02S5000 
-13-2-0-4 i- 3-Q*2r4»^5— 


3d* 9157J JOJ 
3 60 91-69 0 30.1. 
3d« 91360000 
3d*9170i.>724 
3 5o 91140000 

• Jtu. 9^1 3.64 791._ 

38.91260008 
..3d# 912 5d3 9J5 __ 
33*91260000 

Joa.9 L02.362.5_ 

38*91330000 
_3 a 97 4 7 4-l>„ 

36* 913 10000 
3.6.0 90991 5 3 4_ _ 
36* 91310000 
3b* 91.1 7d-b.25._ _. 
36* 912400 00 
30*91351i11 
36*91100000 
_3_£i> 913_?_S-?J3 £l 
38* 91190000 
_ 3d* 9154-^-7-66— 


243.1059000 
..._24 2e9 22 7^-V7 — 
282 * 6374o 00 
_. 282*.7.25.4.0 SO.— 
321*73570 UO 

_^32ilU).594o-63l— 

0 *8454009.0 00 

_1-0^27-bl 6 JU 6.5- . 

3v *9987 Oj 00 

_ 4_0 #.3_4-27_4y 49- ... 

7V * 5 7 16 Ov 0 0 


3J4 o 1 C S3000 

_aa<uu.ii*3S_59- 

273 * 701^000 
__ 2-73u> 7 U2 9-966- 


106*3569000 

^407445^ 


€8o39350GOO 

6 ^3051 


242 *6 17CQOC 
-ZA7: *6227 636- 


32^78560900 


213.53010OC 


182*4299000 
-laXpAASt^AS- 
152 *0246000 
45 :-oQ2 tscax- 


359*3531OOO 


325#» 8239000 
3-Z5A 4X575X5 


118*7110;JO 
J.1Aa4^-UU_0A- 


157 *6965J 00 

_157-4.694U17_ 

197 *4b21^ 00 
ly7 *04b9d 03 
23b *5782J 0 0 
3.5 *3586i 52- . 
2 7 b * p 7 0 1 j ao 
.275* 709-6/19- 


122 ^9580JCC 
_ 9513 - 9 . 11 - 

92*67990000 

-S^Ul729i.7_a4_ 

€2*48060000 
62 o 47 2 5786 1 
22*4041J 0 0 0 

^ix.^A-ooo b oa a -. 


289o6068000 

-23W523A372- 

249*82580041 

_250.^ OA 55A6+2 

205*915200-0 


159*6370000 
160a0415509 
115*4020000 
145^62 71013- 


2o 3 2 3400 00 C 


75*13230000 
-7-5*0-9237304- 


8 


EPOCH (MJDf 


SEMIMAJOR 

AXIS 

(EARTH 

RADII) 

A 


ECCENTRICITY 

E 


INCLINATION 

(DEGREES) 


41 76 4 « -io { 2 9 


4 A 77 1 29 


4 1 77 c * 94i 26 


4 1 79 5 • 5>3 6 2 


4 1 79 2 »96> 24 


4 1 79 ^.94j 9 1 


4 A 60 o«9^j V9 


4 1 61 4*9sj U-J 


4 1 62 0 *94c*5 


4 1 62 7 * 99 ♦ 58 


41 38 4 *97j S»'6 


4 i o4- 1 *94/ 47 


41 *4-6.9** 42 


4 1 65 £ • 9 7 •+ 7 o 


4 1 b5 2,94* 19 


41 S 6 3 *99^ 77 


4 1 87 6 *97^ 60 


4 i. 60 J « 94 / 10 


4 l t>9 0 * _> 29 


4i 69 7 .97 + 29 


41 90 4 47 


4 l 91 1 • 72^36 


1 .1c? 74JC2U 


1«18773cc4j 
1 «. I c 7 y 3 £iii;4_Q__ 
1*1c772cS60 
1*127-720*5.0. 
l-ic77l/wig 
1 . 1 r 77 l Vuli> 

1 • 1 c77C 6 c50 
I.i.l.£77fl.ciii0._ 
l*lE769 91iiO 
1 « 1 6.76J*.*_1,80 

1 a 1 6 7 JC X 6 (j Cl 

_l-*-.l—C-766 » 3 g i ) 

1 a 1 t 768J41 C 

1 « 1^7660490 
— l-»_t c7f.ri0490... 
1 . l€767dC7g 

1 * 1c767 4 76J 

1 « 1 57670800 
...1 ale / 6 60.0_ 

1 a 1 c 7 c 6 564Q 


1 • 167661410 

_l_*.l_c.7.£>6141 0_ 

l aid?657120 
1 * 1 .i 705?12 d_ 
l • 1 6765312 0 


1 # 1 3 7 o 4 9 930 
1*1 S7649930 
laic 7646780 

X. *_1 _£.7 o 46 7. b 0_ 

1 • 16 764 343o 
643430 


1 *1 375369 90 
1*1 c.7 63 g 990 
1*16763j 170 
la 187 o3 317 C 
1 a 1O7o2 7920 
l.*.l.e.7ii32JLi2ii_ 


C* 81662500000-01 3o*9130u000 

J*Jl6o69oOOoD -01 38*91430000 
-0-*. 6 JJj 6_S.43L9.01 D — J 1 3a. 91 4 b 935 >> 
0*01504600000-01 38a91S3u0G0 

- 0*^iXb-.4.2_9-ci96.iiilj^J 1 _< . v) > $ * ? 6 ^ 7 

C ■ 3 1 4 3 *6 OOOOD—01 3d•9165J 00 u 

—C.* 62 J r 66 1 2 6 KP — 0 1 AH . 9 t 1 0 3 ^ ^ 

0 *31450460000-01 38*9179joOO 

.-..C.«-d2_LL27_4 .lAQilQ=ULL_ 3 6* 914542 64 
C*dl35480000D-01 36*51olJ000 

-0*_ttX£A766l 8 DDjrJJLl_3^* 9 1 439663 

C * 2 13 92 70 00 Ob — 0 1 38*91560000 

— 0-«.^.1.2.a_ 9A - f i3J. 3D - U 1 _ . 36 . 91 ct4 > 

0.314 35 7000 00-01 3a*9l530u00 

—0-aJ.O-9.6224QaoO-nl 5a* J1 BKi AS P 

0.3145190000D-01 33*91680000 

—C, .. iljQ 3 95.1 A.ajU-Ql^. 3 o* 92 0 Jl»33 ? 

0.61456500000-01 36*92000000 

—q.*.8.U9-7514 91 40 — 01 38.922^0 074 

0*814 09 oO 000*0 —01 33.92150000 

0*613ol 0 OO 0 OD-o1 36.92390000 

„C .81 5.71d.7l0l_o3.D-01 3 1 ,. 92;? 2 v if ;2 
0*81251700000-01 36.92570000 

_ 0 * d 1 6 2 . 15? fifi l D-Ql .lh^ 42 i A<> 
0*6116000OO 00 — 01 3d* 92 52 o 0 00 
-.0.*.31.76l_Q2.aQ.9D- 0 1 .lri.wyi 7 i ^ 1 

(i • J 1254400000—0 1 38.92570000 

_.a*^JJ^Ilenil7li jr.J-0l 33.9g36a74a 
C.612 A 140000D-01 3o*924*0OOO 




0 *81 27U3O0O00-01 36.92150000 

2_ 3_a ..92 3 7 75.9^. 

0 .d1325cOOOOD-Ol 36.92210000 
--0*0 0 7 1-Z7J3.0.99D -ft 1 . 36 * 925tt 3 3S7 
C.113532000 oO -01 38*92450000 

—£.*606 32 3566 UP-01 3b. 92 752 4*7 
0.61251900000—O1 3d*9246UOOo 
C.«Q_l-OAd45 6 940 -01 3 6 . 92 57 jQ 5 1 
C*6li666000UD-01 36.926I000U 

. _0 . 613 7 76 2 4 0 7D —01 2 6 * 9 2467564 
11.61 0806 GOGOD —0 1 36.92650000 

—0.31 6 l 7 42466D —01 


9 


ARGUMENT OF 
PERIGEE 
(DEGREES) 

CD 


3 13* 142 It# 60 


RIGHT ASCENSION OF 
THE ASCENDING 
NODE (DEGREES) 

a 

a' 

331*9134000 


MEAN ANOMALY 
(DEGREES) 

M 


36.59420006 


354•302oj UO 
334. 7 AHni 9 4 
2 3.445500 OO 

7 3*0 37d O 4 / 00 


3 01 * 62 e 5000 
- 2 -U-L* c 35 d- 5 J_ 3 _ 
271•7429000 


2 4l*3346 COO 


4.868500000 
-4^A51^..S.S9.S9 
3 3 L *52 64-000 
-3 

295^*61 48000 


112.2349JOO 
1 1 2 * 0.6-1 3* 49- 
151* J57brf00 

1 5 a > a j k a o 


190.955 Iv 00 


2 30 * 1 0 9 6 w> 0 0 


<11.2582000 


18L. 1570000 


ISO* 7*73000 


Mpv«.-nnrrii 


120*7056000 


256*6901OOO 


213*44 50000 


rATCT.-VJ 


16 7* 22 96000 


122*5016000 


269 *2015c kjO 
26-9,.*,. 1-5.55/ 21 
306* 70d7j00 


90 *63060000 


60*22390000 


61*521 00600 


44*33760000 


347 »o241c JO 
3_4d.*2„4U.9i 26. 
26 *9 9 220 ^ J 0 


30 *1392000G 


10*38490000 


0*5 49 C 0 U tl C COD- C1 337*9916000 


06 « 60 3o Ui> OO 


229*6541000 


301*8537000 


105*6 57 5j 50 


299•652 60 CO 


263* 5296 0 00 


144 * 88 7 3l/ 0 0 
-L_44..* b 19 3^ 9 a, 
184 • 472jy» 00 


269* 6 1 23000 


239 * 1 1550 C 0 


220*8296000 


174*8277000 


225 *t ?2 1 6 */ 00 
■223*322Zod0 
262*760 IcOO 
09 0^ 43 
JJ 2 » 24 1 4* 00 


209*0453000 

2a9mO_4H0300. 

178*9655000 

17Aa qfc£fih^.S 

148*5643000 


129*7205000 


87*58660000 


to.16460000 


341.277 JJ 00 
d41 * 679 63 45 
20 *42710c OO 




59 *59 33 0^ 0 0 


1 18•4634CC0 


88 *4 0450 00 0 


56 *32790000 


15*96140000 


342* 7159000 


30a*2571 OOG 






SEMIMAJOR 



EPOCH (IVUD) 

AXIS 

(EARTH 

RADII) 

A 

A’ 

ECCENTRICITY 

E 

E' 

INCLINATION 

(DEGREES) 

1 

1' 

4 ■ L 91 8 *97 / 34 

1 • 15 7o£3820 

C »ol034 600 000-01 

3 fa* 9289J0Ou 


1 * 1 C 762 .182 0 

U*dl65l4l5 3j9P_nO 1. 

__3J^*j^2.533 1 03. 

4 i 92 5*$Ji L-6 

1 * lc7ul 6 ji)0 

c«ai ou6 9ooood — oi 

3d.92930000 


1*1 i 7ol^900 

0*6141744954D-01 

3d *_92b9jp.J53.7_ 

4 1 93 2*92-. 02 

1 *1 17o0 j510 

0.61104700000-01 

38.92930000 


1.1 cToC a 51.0._ 

_C_s_6JLl. Z5ZJ. J_3XP™Q.L 

3b. 929 1 JUo? 

4 i 93 9 « 9 7 * 11> 

1 *187O0U430 

0.01133200UOP-01 

3 o. 9287 0 U Ou 


1,lc?6CJ430 

0*o075643136D-0l 

3 b. 96 06 7 1 15 

4 1 ^ 4 * 6 * 9 4 j 0 9 

1*157595150 

C* 611 784O000D-U1 

3b.92760000 


1*127595150 

0.d 0582 C7808D-01 

3 6. 93 10 52_12 

4 1 9^ 3 * 99 j 6 3 

1 . 18756^510 

C . 8 1 £ 26 b 0 OO O D - 0 1 

3b. 92b3oOOO 


_ 1 » 1 d ? 36 5J5LLCL .0 « dPC_7jfl_7_2 .1,6JjjjrPJ_3o,o 913JLXJ-££-£. 


4 1 96 0*97 j 9 6 

41 i-6 7 * 94 * *0 
4 1 97 4 a 3 2-+ 7 5 
4 i 9B 1 a 37 j 14 
4 1 9$ a • 9*+-> 77 
42dO 3 o9jj 66 
4 d 0 i 0 o 7 2d 49 
420 i 7 * S 6 

42 0£ 4 o 93 J \ U 
4203 1 *92j 11 
4203 a *97d 71 
420^5 o $ 5 d c 1 


l»13?b£:541o c. ^ 12o3oouaoO-al 

1 * 1 £75c^4lO- _C 6 €09 35c 1 61 Op--0 1 

1 * 1 c 7 5532 0 0 C.6l03SOC000D-01 

1 , 1 $73652 00 _0* 0117897 5 840 -01 _ 

Uli734l3lO Op €0972600000-01 3bo9297JGQO 

i . J c 754 a 3 1 0 0.dl4691939lD-Ql 30*92682590 

1*137540540 c,d09UOb0d0U0-0l 3b*9£9u00GU 

1 * 167 54o54o - ---- W01-5-309^ 1-0 3 ^-0-1 -38*9253 93^6- 

1*167332040 C*d09299000GD-01 3d*93100000 

1 « 1 6 762 2040 C*~d 1-3912 592 9Q-^0-l— 3d. 9263 34 2-1 — 

1 * 11 752 2 34 0 Codl 07o20000D-01 JO.93360000 

-- 1 -, I ^7-5.2334-0 -— -O*. a-0 7 ..u 9 .q,. 5 , 9 .6W^m — ~^Lfe~-93 

1*127^1^350 0.61ua74GU0oO-01 33*93430000 

UU7 51-936-0 O^^4.94.uy.y5-74i^0l —5 a * 9 3773 4 31 —- 

1 o 1c?512350 C * d 1 1 095 uC u JU-0 I -Jo.93520000 

i *ld75L235u — -C *-005515431-VU-—0-1-3-^^93-04-50-7-0— 

1.1<373Cc 330 0* 21 1 01 00 0000-01 3c‘. 935/0000 

—1-*-4 6 7 454^3-30-C-W03 24d 3 ,1 ,J 2D^4l-33*-93 72 0 L 4 o 

1.15 7 SO 1 Cb 0 C*610l^ZOOOOD-Ol 3b.*3510000 

1 * -L-S 75 01 C 6G 0*diT S3.Q31ji_lJ.r--iLl-lt4-934-345it 4 . 

I*16?4v20i0 (,* 0C662 dUOOOD-O 1 36.93310000 

1 -*-l ~ 74 9 J o 10-- - - C »-0-I3o6 7_4_d22.O-*-‘0Ll-3i3-o- * 3 Q2 i 55 2 

1« l€74o4 740 C .0O96b4oGO0U-01 36.93340000 

... .l+}.*r3 4 feu 7 4 0 0 O 3 1 39 1 j 3 7G 4 2-U J-3d* 9 2 97 0 52 3 


36.92730000 
_3_fl_* 92 jS0_+j>_9 0 
36.9267O0U0 
36 * *2 7 3 T 0 39 


440-5 2 o iii j 96 
40 05 9*974 U 


1.1c74730 3 o C* 80967200009-01 3t.93300000 

1 * 16 7-4 7^u.9-0_C-*31-4 j55<LJ4 040.t-0 1 3 6. 93 02 9 2 74 . 

1.127472960 C.61 0 1 U2QOO-OU-U1 36*93460000 

1*1 6 74 72 9b(J - G. 6 Li -06 34 73-uD“G-l _ 3o. 93 4 J.43*-6 


ARGUMENT OF 
PERIGEE 
(DEGREES) 
co 
u) r 

99 * 1 629 Jw- 00 

_ 99 .0 89 1-9.6.67. . 

13d *34030 00 
13d*0 0l3i 6.8 
17 7.6243^00 

_*77 pQ690rf ofl_ 

217*006 1 j 00 

_ _2_l_b_*_ 7 535^5 3_ 

256*2033^ DO 

2 5 6.1 Oo 4 * 36 
29 5 *4 76 £u GO 

33 4.623 To 00 
335 * 2051966 
1 3 .9200 Oj00 
_1 4.346 45^ 46__ 

53.1 0 190j QO 
S3o 37677i: £0 
92.72150^00 

-02-^69+47-=. 2 7- 

131o6876j 00 

-1-3-4—0-5-02 — 

216*17530 00 

- 2 1 5 *8 360035 

255*27570 GO 

-.,.^5^-1- 722^96 - 

2*4*767^300 

333.8350300 

- 3 3 4 * 2 0 3 73- 7-3-- 

13.0 1 850l/ 00 

52 *6^3 50j 0 0 

- 52-0- 9 / 6 2 / ■ .> 6 A- 

V1*665203 00 

- 9 1 . 6 5 0 l 0 / 4-3— 

130 «9 603J00 

-130* 6 610239 . 

170*5143300 


RIGHT ASCENSION OF 
THE ASCENDING 
NODE (DEGREES) 

a 

a’ 

27*92700 OC C 

357.6555000 

_3 E7. _EAS-9.544_ 

327WSC4QC0 

227 o- 7L6.30.521_ 

2 97* 39 60 0-00 
£97 * 39C1945_ 
267*3218000 
_2_6_7_-_32 G.0 Z40.„ 
236*9111000 

_ Z2J^9±AJl±^— 

206*6260000 

176*760000-0 
176^ 7 61 1.5 6 0 _ 
146*6732000 
146* 677 6t E3 
116*2^93000 

-1- 1 1 ^£-* 8 4 4 8-9- 

66*19880000 

£1*46030 0140 

2 1 *4 7 4 » 3 g . 4 S 

3 5 1 *4 07 1000 
- 36 i4>4 05262 4- 

321*0003000 

- £2-1-* 00333^9-' - 

290*9179000 
89 0 *9 34 4 496 
260*627 1 000 

- 240^ 6 4 4 2 794— 

£30*4327000 

- 200-*-437^W5— 

200*3521000 

- jflO *‘ 3 S c 5 571 — 

170*2741GCC 

129*6712000 
.-12 9 * 6 - 6 - 39399 — 


MEAN ANOMALY 
(DEGREES) 

M 

M ; 

270* 1 737-000 

.- 270^24 7fe265v. 

228*2052000 

__220*5452343 

1 82 ■* 9772000 

-__1 83*.A t 36.5i3ft 

137 * 0063 000 


137 ft A2Qflafe2 


94^69310000 

56^39830000 

56*22216035 


2 t«46650000 
21o10564031 


346x. 26 03000 
3 4 7 n fl3 043 (H- 


314m1335000 
3134 35 77847 
276*6451000 
276*6671933 


235*3905000 

^36«6962222- 

138*1221000 

46 04 59*^- 


$5*66000000 
7 6 3 6 3 3 4 7 


57*0465000^0 


22*30000000 

-a- i - ^ o i -fo p fraa n 


369*0266000 
-36 0 -^ 0 $- 65& S3- 
314*4670000 
-34-6o 20 6 3 - 7 -60 ~ 
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Figure 2. Resonant Variation of the Inclination for the Orbit of 1967-14F 



The residuals in inclination (rms) from this trajectory (and others including 
those with resonant effects) are shown in Table 2. The calculated trajectories 
were "fit" to the mean element "observations" by differential correction of 
initial elements and other model parameters (see i.e. , Morrison, 1970 or 
Wagner and Douglas, 1970) under the conditions stated in Table 2. The principal 
difficulty in these orbit determinations was in following the mean anomaly of the 
satellite which (with a perigee of 580km) underwent fluctuations, of tens of de¬ 
grees from drag error. The level of this error was less than 10% after correc¬ 
tion for a single drag coefficient. Yet it was necessary to include additional 
accelerations in the semimajor axis and mean anomaly to permit the critical 
resonant longitude be calculated to better than 5° (rms), (This error is 

consistent with the formal accuracy of the resonance determination.) 

A more readily available data source for satellite objects, the North Ameri¬ 
can Air Defense Command’s SPADATS elements were also evaluated for 1967- 
14F by the same mean element program. (See run 2 of Table 2,) Figure 3 
shows the inclination residuals from a nonresonant trajectory with this data over 
a somewhat shorter arc in 1972-1973. For most of the period no detail at all 
can be seen. At the end of the arc some definition of the resonance appears. 

The NORAD-SPADATS data quality is clearly not uniform in this arc. The 
small "acceptable" portion could not be used for adequate resonance recovery. 
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Table 2 


Results of Orbit and Coefficient Determinations for 1967-14F Using 

Navy and NORAD Data 


Run 

Field Used 

Residuals * 
in 

Inclination 

(rms) 

<10- 3 deg’s) 

Data Used 

Comments 

i 

SAO SE 2 (non 
resonant) 

5.41 

NAVSPASUR - 
All elements 
MJD 41303 - 
41989 

Drag, radiation, a 
and M coefficients 
solved from data. 

C ommensurability 
at 41985. Reso- 

inance i i"I"clear, 

1 

residuals about 
0 .001°. 

2 

SAO SE 2 (non 
resonant) 

6.25 

SPADATS- 
All elements 

MJD 41302 - 

41936 

Same as above. 

Resonance in I not 
seen except pos¬ 
sibly after MJD 
41700. Residuals 
about 0. 01°. 

3 

SAO SE 2 (non 
resonant) 

6.83 

Same as 
run #1, but to 
MJD 42060 

Same as run #1, 

4 

SAO SE 2 + 3,0 and 
(23, 13) solved from 
data 

109(C,S) 23>13 - 
(-20.7 ± 1.4, -69.8 
± 1.2) 

1.13 

Same as 
run #1. 

Same as nm #1. 
Weight on I: 

0.0002°. Maxi¬ 
mum along track 
error: 6°. 


*0bserved - calculated values from converged mean element trajectory 
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Table 2 (continued) 


Run 

Field Used 

Residuals * 
in 

Inclination 

(rms) 

<10 -3 deg's) 

Data Used 

Comments 

5 

SAO SE 2 + 3,0 and 
(23, 13) solved from 
data 

10 9 (C, 8)23,13 = 
(-22.0 ±1.2, -71.8 
± 1.1) 

1.29 

NAVSPASUR - 
All elements 
to MJD 42060 

Drag, radiation, a, 
a, MUl, M^i , 

MI 5 ! coefficients 
solved from data. 
Short period lunar 
terms used. Max¬ 
imum along track 
error: 10°. Cor¬ 
relation coefficient: 
(C,S) 23;13 = -0.38. 
Weight on I: 

0 .0002°. 

6 

SAO SE 2 + 3,0 and 
(23, 13) solved from 
data 

10 9 (C,S)23,13 = 
(-21.2 ± 3.4, -71.3 
±2,4) 

— 

1.29 

NAVSPASUR - 
Inclination 
data only - 
to MJD 42060 

All state elements 
except "I" and ra¬ 
diation coefficient 
fixed from solution 
in run #5. Corre¬ 
lation coefficient 
<C,S) 23) i3 = 0. 34. 


*Observed - calculated values from converged mean element trajectory 
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Figure 3. Inclination Residuals from a Nonresonant Trajectory for 
1967-14F Using NORAD-SPADATS Observations 




RESONANCE RECOVERY 


The geopotential indices for this (mean) primary resonance (C-2p + q=l) 
are m = 13, q = 0 and all £ > 13 for which 

£ = 2 p + 1, p < £. 

Thus the £ are all odd and any one of the (£, m) terms (13,13), (15,13), (17,13) 

. . . can be recovered from the resonance in inclination on this orbit. 

It is noted that the first of the side band primary resonances (q = ±1) which 
are affected by the even degree geopotential term (and m = 13) have only slight 
effect in this data span since to = 5.6° /day (relatively large compared to i// )3i0 
in this period). The q = -1 resonance has a period of 180 days at 4300MJD de¬ 
clining to 60 days at 42060 MJD. The q = +1 resonance has a period of only 40 
days at the beginning of the data span, increasing to 70 days at the end. In ad¬ 
dition, the amplitudes of these terms are of order e' q , where e is the orbit's 
eccentricity, [see Kaula, 1966, p. 37J or reduced by a factor of 0. 08 with re¬ 
spect to the q = 0 resonance for 1967-14F. 

The fact that (essentially) only a single geopotential term can be recovered 
from each resonance is a consequence of the fact that only a single harmonic 
perturbation is responsible for it (see King-Hele, 1973a or Wagner, 1973). The 
amplitude of this term changes only as a consequence of the change in frequency 
through the resonance. The fundamental amplitude (a weighted sum of the res¬ 
onant geopotential coefficients) remains constant. The scale of this sum is 
determined by the actual resonance perturbation. However, even for the same 
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fundamental frequency ( \jj ) this weighted sum is different for 4 of the Kepler 
element variations (a, e and "I" have the same sum, but co, £1 and M are 
distinct). 

Why (physically) only 4 of the 6 sums (or amplitudes) yield independent in¬ 
formation is unclear. It may be due to the choice of the classical Kepler ele¬ 
ments to express the perturbations. In particular, the choice of the mean 
anomaly instead of the true anomaly introduces an infinite number of frequencies 
to express the effects of a single geopotential harmonic on an orbit. In some 
sense this must dilute the information content of any single frequency. But in 
any case, for the dragged orbits, only the inclination variation (essentially free 
of drag error) has provided unambiguous recovery of resonance information. 
However, recently, King-Hele (1973c) and Winterbottom and King-Hele (1974) 
have shown that near circular orbits can suffer significant resonances of £2, e 
and co apparently quite distinct from drag (and other) effects. A strong reso¬ 
nance in e was also seen by Wagner (1973) in the very slow decay of the Vanguard 
3 orbit. On the other hand no significant resonance of other elements on the 
DIADEME 2 fragment has yet been seen, though the full passage will not be over 
till late 1974. 

I chose the geopotential harmonic (23,13) to absorb the inclination reso¬ 
nance on 1967-14F and at the same time made a 5% adjustment of (3,0) to cor¬ 
rect odd zonal model error (see Table 2, runs 4-6). The same differential 
correction program was used in these adjustments as previously to reveal the 
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resonance. But added weight was given to the inclination data in these adjust¬ 
ments and additional secular accelerations were used to reduce the along track 
or phase error of the resonance. Full data correction runs were made for 686 
and 757 day arcs (runs 4 and 5). The inclination residuals were significantly 
larger in the longer arc, possibly showing the influence of the secondary reso¬ 
nance with m = 26. But this conclusion must wait till the passage is complete 
and the smaller (higher frequency) effect can be well separated. In any case, 
the (23,13) recovery is not substantially altered over the longer span. A final 
run was made using the inclination data only, to confirm that the other elements, 
influenced by drag, were not distorting this result. The (23,13) recovery in 
this correction (run 6) was between the values in the two full data analyses. 

The solid curve in Figure 2 shows the computed resonance in inclination 
from the best (most representative) results to date (run 5). The values for the 
harmonic coefficients themselves are somewhat larger than Kaula’s rule [Kaula, 
1966, p. 98j . But this is not significant since it is only a linear sum of reso¬ 
nant terms which is well determined. 

RESONANT CONSTRAINT (LUMPED COEFFICIENTS) 

Following the method developed by Gooding (1971) and elaborated by Wagner 
(1973) the resonant inclination variation for 1967-14F (a = 1. 88 e.r., e = 0.082, 

I = 38.92°) is determined by the lumped sine and cosine terms: 
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(C,S) 13 ,o = 0.023(C,S) 13 B -0.172{C,S) 15 13 +0. 505(C, S) 17 13 
-0.884(C,S) 19>i3 + 1.000(C,S)2i,i3 -0. 673(C, S) 23 ,i3 +0. 099(C, 8)25,13 
+ 0 . 295(C,S) 27 , 1 3 -0.279(C,S)29,13 + 0 . 018(C, S) 31j13 +0.156(C, S) 33 ,i 3 

-0.105(C,S) 35 13 -0. 036(C,S) ?7 l3 

+0.085(0, S) 39jJ3 -0.021(0, S) 41>13 + . . . (1) 

The terms in this series are for fully normalized geopotential harmonics 
[Kaula, 1966, p. 7). The weights are just the fundamental amplitudes (without 
the rate denominator) of the linear perturbation of the inclination due to a fully 
normalized harmonic [Kaula, 1966, p. 40J. The sum of these fundamental 
amplitudes (with C and S coefficients) are merely the coefficients of the cos 
(i// )3 0 ) and sin (>p yi 0 ) terms determining the rate of the resonance variation of 
the inclination. It is these two (lumped) coefficients which are actually "ob¬ 
served 1 ' in this resonance. The reason (23,13) was chosen to absorb the effect 
was because it had a high weight in (1). It also was the lowest degree 13th order 
term not present in the SAO SE 2 field. Using the (23,13) weight of -0.673, the 
"observed" lumped coefficients are (best results): 

10 9 (C,S) 130 = (14.8 +0. 8, 48.3 +0.7), (2) 

with a correlation coefficient of -0. 38. 

With regard to "lumped coefficients", if their perturbations can be properly 
identified, the weights of their constituents can serve to extend the information 
in them to any degree. Such "lumped coefficients" have been reported many 
times from "shallow resonant" satellite orbit analyses [i.e., Yionoulis, 1965; 
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Murphy and Cole 1968; Gaposchkin and Veis 1967; Douglas and Marsh, 1970]. 
One difficulty with ’’shallow resonant coefficients” is that except for near circu¬ 
lar orbits, the side band resonances (i. e. , q = ±l)will also be observable. They 
will not be separable unless long arcs of data are analyzed spanning at least a 
rotation of perigee. But if the "lumped coefficients" are well determined (and 
identified), as here and in the previous 15th and 11th order analyses, they can 
serve as absolute benchmarks for geopotential determinations. An example of 
this use has been given by Wagner (1973) for a previously poorly observed 11th 
order resonance. Here, I calculate the "lumped coefficient" for the 1967-14F 
resonance from Equation (1) with a number of recent fields which are well repre¬ 
sented by 13th order terms. The results are presented in Figure 4. The merit 
of each of these fields is now immediately apparent when compared to the "ob¬ 
served" coefficients. 

DISCUSSION OF RESULTS 

The "point" in Figure 4 representing the "observed" lumped coefficients 
has two uncertainty circles about it. The inner one represents the formal 1 a 
uncertainty of the "best" solution. The small solution shift from the abbreviated 
data analysis is also seen. The larger circle represents the uncertainty in the 
lumped coefficient if terms of degree higher than 30 are ignored (as they gen¬ 
erally are in current global solutions). This expected "truncation error" is 
calculated as the root sum of squares for the neglected terms in Equation (1) 
using Kaula’s rule (10~ 5 /£ 2 ) for the harmonic coefficients (Table 3). It is also 
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CALCULATED FROM 
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Figure 4. Lumped Harmonic for 1967-14F Resonance 



Table 3 


Estimated Cumulative Effect on Lumped Harmonic for 1967-14F, 

From Geopotential Terms 



8 

13 

15 

17 

19 

21 

23 

25 

27 

_ 

29 

31 

33 

35 

37 

39 

41 

10 9 RSS* (all terms > C) 


40.9 

40. 9 

40.2 

36.2 

26.6 

13.9 

_ 

5.8 

5, 5 

3.8 

cn 

* 

1.9 

1.3 

r - —i 

0.9 

0.8 

0.6 


*Using Kaula’s rule and Equation (1). 



of interest to calculate the a priori variance of the total lumped coefficient, as 
its square root can be thought of as a reasonable (la) value for it. Using Kaula’s 
rule again [ over all the coefficients in Equation (1)] this value is indicated by 
the circular arc centered on the origin in Figure 4. 

The convergence of the series in Equation (1) is actually stronger than it 
appears. The weights can be shown to decrease (on average) at least as fast as 
l/G while the coefficients behave as 10" 5 /2 2 . Comparison of the observed 
lumped harmonic with its a priori standard deviation shows its value is slightly 
greater than Kaula's rule would have predicted. In terms of Kaula’s rule, the 
formal statistics say the lumped harmonic is determined to better than 1 part 
in 40. 

The fields chosen for comparison with the ’’observed" harmonic all used 
basically the same shallow 13 th order resonant orbits but with different amounts 
and kinds of tracking data. They all extend to at least the 22nd degree in 13th 
order terms. The purely satellite data fields are GEM 5 and PGS62 (F. Lerch, 
Personnal Communication, 1973 and 1974). The others are combination solutions 
with surface gravimetry data. The Smithsonian fields (GaposchkinandLambeck, 
1971; Gaposchkin, 1973) used analytic techniques for both orbit and geopotential 
determination. The satellite data in these consisted only of camera and laser 
observations. The fields originating at Goddard Space Flight Center employed 
numerical integration to rationalize the orbital data. GEM 5, 6 and PGS-62 con¬ 
tain significant amounts of electronic data (radar range and range rate, and 
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Doppler observations) on the 13th order shallow resonant orbits. PGS 11 (F. 
Lerch, Private Communication, 1974) contains only camera data for the satel¬ 
lite observations. 

As for the truncation of 13th order terms, GEM 5, 6 and SAO SE 2 extend 
to (22,13), SAO SE 3 goes to (23,13), and PGS 11 and 62 include all 13th order 
terms to (29,13). 

Yet in spite of all these differences, the fields are clustered fairly closely 
in Figure 4, showing the dominance of the similar" 13th order satellite 
information. The Goddard fields (with substantially more data) are somewhat 
closer to the observed harmonic, but the truncation is significantly different for 
them. In fact the (probable) truncation error alone would account for all of the 
distance of the GEM 5 and 6 models from the observation. On the other hand 
the SAO SE 2 (at the same truncation) is significantly farther from the observed 
value. (It is somewhat beyond the probable truncation error for all the terms of 
degree greater than 22.) The SAO SE 3 (with 23,13 terms) is marginally closer 
than SAO SE 2 to the observation but the added terms should have made a much 
greater improvement still. The same situation holds (and even more strongly) 
for the more recent Goddard fields (PGS 11 and 62) which extend to (29,13). 
Here, no improvement is seen over the earlier GEM solutions but the truncation 
error (given by the dotted circle about the observation) should be much reduced. 

The simplest way to interpret these comparisons is to say that the higher 
degree terms (i.e. , above about 21) for this resonance are not well determined. 
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The conclusion then follows that these terms will be significantly improved with 
the use of the lumped harmonic (i. e., Equation (2)) for the DIADEME 2 frag- 
ment. The simplest example of such use would be to add a single high degree 
term to an existing field. Of course this would not be a realistic solution on 
two counts. First, it would upset the 13th order ties (correlation) in the existing 
field. Secondly, it would ignore the contributions of still higher degree terms. 
But it does give a representative value (exact for 1967-14F) and a first approxi¬ 
mation of the realistic numbers. Using GEM 6, I find the added (23,13) term 
(C.S) to be 10“ 9 (-16. 9, 0.3) which has an ms of 12. 0 x 10' 9 compared to 18. 9 
x 10' 9 for Kaula’s rule. 

SUMMARY AND CONCLUSIONS 

A strong 13th order resonance has been observed and analyzed from U.S. 
Navy Tracking Data on the slowly decaying orbit of a DIADEME 2 fragment 
(1967-14F). The exact commensurability for the orbit occurred in late 1973 
and the major changes due to the resonance will be over by late 1974. Never¬ 
theless, apparently stable and well determined values of a lumped harmonic 
for this resonance have been found which should significantly improve 13th order 
geopotential terms to about as high as degree 33. There is fairly close agree¬ 
ment of calculated values from recent fields with the '’harmonic" observation 
(within 20%). The major part of the discrepancy is probably due to poorly 
known coefficients above degree 21 in these fields. 

Substantial improvement of 13th order and high degree terms will be seen 
with use of the lumped values (and the linear constraint) in combination solutions 


with other data. 
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